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O . Abstract 

^ \ Semiclassical analysis of the shell structure for a reflection-asymmetric deformed os- 

cillator potential with irrational frequency ratio uj^_/uJz = V^ is presented. Strong shell 
O ' effects associated with bifurcations of short periodic orbits are found, which occur for a 

^ • combination of quadrupole and octupole deformations. 

> 

X 



t Submitted to Physis Letters B. 

■''E-mail: aritaOruby . scphys .kyoto-u. ac . jp 



Study of reflection-asymmetric deformation is one of the exciting current subjects in 
both nuclear structure and micro-cluster physics [1-4]. As is well known, nuclear deforma- 
tion is intimately related with the single-particle level structure. The system favors such 
shapes that exhibit strong shell effects and make the level density at the Fermi surface 
lower. We have analyzed in Refs. [5-7] shell structures in a superdeformed (SD) oscillator 
potential with the octupole deformation added. The model Hamiltonian adopted is 

H=^ + T. ^^ - XsoMu;',[r^Y,om", (1) 

i 

where the double primes denote that the variables in the square bracket are deflned in 
terms of the doubly-stretched coordinate x" = {ijJi/ujQ)xi. The frequency ratio uj^/oJz is 
taken to be 2 with uj^ = ojy = uo^^. At A30 = 0, the shell structure for this Hamiltonian 
is mainly determined by the flgure-8 type family of periodic orbits. However, the planar 
orbit family in the x-y plane having shorter period give rise to an interference effect on 
the shell oscillations and result in modulation pattern in the smoothed level density, called 
supershell. We found that the supershell effect signiflcantly develops at octupole deforma- 
tion parameter A30 — 0.4 [^. This is explained as due to the changes in stabihties of the 
two classes of periodic orbits against the octupole deformation |^, |^. 

As a continuation of the above work, we have carried out semiclassical analysis of the 
shell structure for the Hamiltonian (^ varying both A30 and the axis ratio ujj_/uJz. The 
details will be reported elsewhere ||^. As an illustrative example, we discuss here the 
case of an irrational axis ratio oJ^joJz = "\/3- Figure |l| shows the single-particle spectrum 
calculated as a function of A30. There is no prominent shell structure at A30 = because 
of irrationality of the frequency ratio. However, a signiflcant shell structure appears at 
A30 — 0.3. The oscillating level density smoothed by the Strutinsky method is shown in 
Fig. ^ One clearly sees a prominent shell oscillation. Evidently, this shell structure arises 
as a consequence of the combination of quadrupole and octupole deformations. One may 
also notice a supershell structure, which is found to be associated with the interference 
between classical periodic orbits with the periods T ~ 27r/uj± and 2tt/uJz. 

In classical mechanics, there are only two types of periodic orbits at A30 = 0; one is the 
Isolated Linear orbit along the z-axis (IL) and the other is the continuous family of elliptic 
orbits in the x-y plane. Among the latter family, only the Planar-A type (PA) orbits (in 
the plane containing the symmetry axis) and an isolated circular orbit survive when the 



octupole deformation is added. In Fig. ^, we show some short periodic orbits for A30 = 0.3. 
It is noteworthy that orbit PB appears at A30 ^0.29 due to an isochronous bifurcation of 
orbit IL. Likewise, orbits PC and PD are born by a saddle-node bifurcation at A30 — 0.28. 
Orbits PE and PF arise from a period-tripling bifurcation of orbit PA at A30 — 0.2. To 
show how these bifurcations occur, we plot in Fig. ^ traces of the monodromy matrices 
for these orbits as functions of A30. The monodromy matrix M describing linearized 
dynamics about the periodic orbit is defined by 



^'■^^^°^- dZ^{t = Q) 



(2) 



where Z denotes the phase space vector (p,q), Zr^■, a point on the periodic orbit r, and 
Tr the period of the orbit. Z^{t) represents the components orthogonal to a manifold 
formed by the continuous set of the periodic orbit, and here has dimension 2 (or 4 for an 
isolated orbit). The eigenvalues of M are independent of Zr^ and, due to the symplectic 
property, they are represented as (+/— )(e°', e~"), a being real or pure imaginary. Stable 
orbits have imaginary a and |TrM| < 2, while unstable orbits have real a and |TrM| > 2. 
The periodic orbit bifurcation occurs when the eigenvalues of M become unity, namely, 
when Tr M = 2. The bifurcations mentioned above are clearly seen in Fig. ^. 

In order to see the feature of the bifurcations in detail, let us examine the Poincare 
surfaces of section. Since the system under consideration has axial sjTumetry, it reduces 
to a two-dimensional system with cylindrical coordinates (p, z) having a fixed angular 
momentum p^. Figure ^ shows the Poincare surfaces of section {p,Pp) defined by 2; = 
and pz < 0, for p^ = 0. The orbit IL (corresponding to the origin) is stable at A30 = 0.28 
and accompanies tori about it. These tori are significantly distorted until a saddle-node 
bifurcation occurs at A30 — 0.283. In the figure for A30 = 0.29 we thus find a new phase- 
space structure that accommodates two emergent periodic orbits; the stable orbit PD 
(corresponding to the pair of islands) and the unstable orbit PC (the pair of saddles). At 
A30 — 0.292, an isochronous bifurcation of orbit IL occurs, generating the stable orbit PB 
(corresponding to the new pair of islands seen in the figure for A30 = 0.3), and the central 
torus becomes unstable. 

Now let us recall the correspondence between quantum and classical mechanics. In the 
Gutzwiller trace formula, the density of levels g{E) = J^n^i^ ~ ^n) is represented as a 



sum of contributions from classical periodic orbits [10 



g{E) = -g{E) + Y. Anr{E) cos (^^M^ _ 1^\ . (3) 

The first term on the right-hand side is the Weyl term which is a smooth function of 
energy. The second term is composed of the periodic orbit sum (r denote each primitive 
periodic orbit and ^(7^0) the number of repetition of the orbit) and represents the quantum 
correction to the first term. S'r = f^p ■ dq is the action integral along the orbit, and /i„r 
denotes the Maslov phase. Our aim here is to analyze gross features of the single-particle 
spectrum, namely, the shell structure. Let 5E denotes the energy resolution which we are 
interested in. Then the periods T = dS/dE of contributing periodic orbits are restricted 
by the following relation: 

A^ = S{E + 6E) - S{E) ^T-6E< 27r^, (4) 

namely, by 

Thus, we need only to pick up these finite number of short periodic orbits. Next, let us 
discuss the amplitude factor Ar{E) in the trace formula. For systems in which chaos is 
fully developed, periodic orbits show strong instability and are well isolated in the phase 
space, so that the stationary phase approximation (SPA) seems to work well. Then the 
trace integral can be performed by the SPA with the result 

Anr{E) = -\ ^'- =, (6) 

^^Y/|det(l-M/)| 

where Mr denotes the monodromy matrix and T^ the period of the primitive orbit r. This 
expression may be good when all eigenvalues of M^" are sufficiently far from unity. At 
the point where a pair of the eigenvalues becomes 1 (a bifurcation point of the orbit), 
this amplitude factor suffers divergence. It is due to the break down of the SPA: In the 
coordinate expansion about the stationary point, one of the coefficients of the quadratic 
terms vanishes and higher order nonlinear terms become important. Such a situation 
occurs quite often in mixed systems, where one observe admixture of regularity and chaos. 
This is an unsolved difficult problem of quantizing soft chaos. 



To find the link between the quantum shell structure seen in Figs. |I|, ^ and the properties 
of the periodic orbits, let us consider the Fourier transform of the single-particle level 
density [0 



oo r^isE 



F{s)=l dE^giE). (7) 

"'0 \/E 



To remove the ambiguity associated with the cutoff in energy, we adopt a gaussian convo- 
lution with /(x) = exp(— x^/2) as 



FAs{s) = Jds'f{{s-s')/As)F{s'). 



Using the scaling property of the system under consideration, H{ap,aq) = a'^H{p,q), 
the quantum and the semiclassical expressions for Fas(s) become 

FLT\s) = E4^/(^n/i5max), (9) 

n V -^n 

F^hs) = FUs)+Y.^nr{l)e'^'--"f{{s-nTr)/M. (10) 

nr 

where the cutoff energy is defined by -Emax = h/As. The semiclassical F^'^^\s) has the 
functional form with successive peaks at the periods of classical periodic orbits, and the 
height of each peak represents the intensity of the corresponding orbit. Thus, one can 
extract information about classical mechanics by calculating the quantity (P) and com- 
paring it with the expression ([T0|). One may expect that the singularities in the classical 
phase space due to the periodic orbit bifurcations will affect the properties of the quantum 
system in a critical manner. 

Figure ^ shows absolute values of the Fourier transform (^ as a function of both s 
and A30. It is seen that the peak at s ~ ^/3 (in unit of T± = 27r/c<j^) significantly grows 
up as A30 increases and reaches the maximum at A30 = 0.34. At small A30, this peak is 
determined by the linearized dynamics about the orbit IL. As A30 increases and approaches 
the bifurcation point A30 — 0.3, the peak becomes higher as expected from Eq. (^. In 
addition, the resonant torus (which brings about orbits PC and PD) also contributes to this 
peak because of the similarity of the frequencies. Furthermore, three-dimensional periodic 
orbits (having periods similar to them, not shown in Fig. ^ contribute to it too. Thus 
these bifurcations of orbits that take places at almost the same value of A30 are responsible 
for the strong enhancement of the Fourier peak. 



One should also note that the peak at s ^ v3 takes the maximum, delaying somewhat 
from the bifurcation points A30 = 0.28 ~ 0.3. Another peak at s ~ 3, corresponding 
to orbit 3PA (triple traversals of orbit PA), reaches the maximum at A30 — 0.28, also 
showing a delay from the period-tripling bifurcation point A30 — 0.2 of orbit PA. Similar 
delay phenomena were noticed in our previous analysis for the SD case with the frequency 
ratio Lo^jiOz = 2 [^]. Thus, we have found a presumably general tendency of the Fourier 
amplitude near bifurcations: The amplitudes exhibits significant enhancements due to 
the periodic orbit bifurcations, with the maximum point being somewhat shifted into the 
postbifurcation region. It would be very interesting to investigate these behaviors of the 
Fourier amplitudes in more detail. 

In conclusion, we have made a semiclassical analysis of the shell structure for a 
reflection-asymmetric deformed oscillator. Prominent shell effects arise for a combination 
of the quadrupole and the octupole deformations. It is shown that these shell structure 
can be understood in terms of short periodic orbits and their bifurcations. It would be an 
interesting future subject to investigate how the new shell structure discussed in this pa- 
per persists in other reflection-asymmetric deformed potentials having radial dependence 
different from the oscillator- type. The author thanks Professor Matsuyanagi for many 
valuable discussions and carefully reading the manuscript. 
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Figure 1: Single-particle spectrum of t he Hamiltonian (|ID with the frequency ratio 
^±/^z = v^, plotted as a function of A30. Dashed and solid lines represent levels with 
zero and nonzero f^-quantum numbers, respectively. Each solid line is composed of two 
degenerate levels with K and —K. Magic numbers appearing at A30 ~ 0.3 are indicated 
taking the spin degeneracy factor 2 into account. 



Figure 2: Oscillating level density gosc{E) = g{E) — g{E) for A30 = 0.3. The Strutinsky 
energy smoothing is done to a width 5E = huj±/2. The numerals indicate the magic 
numbers corresponding to the shell gaps. 



Figure 3: Some classical periodic orbits for the Hamiltonian ([T|) with the frequency ratio 
'-^±/^z = a/3 and A30 = 0.3. 



Figure 4: Traces of the monodromy matrices in the neighborhood of the bifurcations 
of orbits IL (upper panel) and of orbit 3PA (lower panel). The bifurcations occur when 
TrM = 2. 



Figure 5: Poincare surfaces of section {p,Pp) defined by 2; = and pz < 0, for p<p = and 
A30 = 0.28 ~ 0.31. 



Figure 6: Fourier transforms (|^) of the quantum level densities with -E'max = lO^c^o, 
plotted as functions of s (in unit of 27[/u!±) and A30. 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9407004vl 
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This figure "fig2-l.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9407004vl 
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This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9407004vl 
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This figure "figl-4.png" is available in "png" format from: 
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